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Abstract. We analyze the dispersive properties of a Dirac system perturbed 
with a magnetic field. We prove a general virial identity; as applications, 
we obtain smoothing and endpoint Strichartz estimates which are optimal 
from the decay point of view. We also prove a Hardy-type inequality for the 
perturbed Dirac operator. 



1. Introduction 

The Dirac equation is one of the fundamental systems of modern physics and 
mathematics, used to describe a spin 1/2 particle in quantum electrodynamics. 
Even if the physical interpretation of Dirac fields is not completely unambiguous, 
the very rich mathematical structure of this system makes it an interesting object 
of study. We refer to [TH] for a thorough treatment of the subject, including the 
physical validity of the model. 

We fix our notations. The Dirac equation is the 4x4 constant coefficient system 



(1.1) 

where it : R t x 



iut = m(3u + T>u, m G I 
Z 4 , the operator T> is defined as 

3 



fc=i 

and the 4x4 Dirac matrices can be written 



OLk 





O'k 



<?k 





0- 



in terms of the Pauli matrices 
h = 



1 
1 



02 = 



1,2,3 







CT3 = 



1 



The coefficient m is called the mass and we shall distingush the massless case 
m = from the massive case m ^ 0. This distinction is important in relation to 
the dispersive properties of the equation. Indeed, the commutation rules 



aiak + akat = 25 k il 4 



imply the identity 

Thus we have the property 

(idt 



V 1 = -AJ 4 



V)(id t +V) 



(A - d\ t )U 
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showing the intimate relation between the Dirac and the wave equation. A similar 
computation in the massive case produces a Klein-Gordon equation with positive 
mass to 2 : 

(id t -V- m(3)(id t + V + m(3) = (A - m 2 - d 2 t )h 

It is then straightforward to derive dispersive and smoothing properties of the free 
flows from the corresponding ones for the scalar equations using these identities 
(see e.g. [5|). Our purpose here is to extend these properties to the case of system 
perturbed by a magnetostatic potential 

A = A{x) = (A 1 (x),A 2 (x),A 3 (x)) : R 3 — > R 3 . 

In virtue of the principle of minimal electromagnetic coupling, the influence of the 
field is introduced in the equation by replacing the standard derivatives with the 
covariant derivatives 

V A ■= V - iA, 
so that the operator T> is replaced by 

3 

(1.2) V A = i- 1 Y,Mdk~iA k ) 

We shall also use the unified notation 

(1.3) H = r 1 a- V A + m(3 = V A + m(i 

for the perturbed operator, which covers both the massive and the massless case. 
Thus our main goal here is to investigate the dispersive properties of the flow 
u = e ltH f relative to the Cauchy problem 

(1.4) iu t (t,x) +Hu(t,x) = 0, u(0,x)=f(x) 

where u(t, x) : R x R 3 -> C 4 , f(x) : R 3 — > C 4 , and to e R. 

It is natural to require that the operator TL be selfadjoint. Several sufficient 
conditions on the potential A are known; e.g., if the field A is smooth or satisfies 

\A{x)\ <A- + 6, a < 1, b> 

then Tt admits a unique selfadjoint extension. We refer to [19] for a discussion of 
this problem; here we prefer to make an all-encompassing abstract assumption on 
the operator: 

Self-ADJOINTNESS ASSUMPTION (A): the operator H is essentially 
selfadjoint on C C °°(R™), and in addition, for initial data in C C °°(R"), 
the flow e ltn f belongs to C(R, H 3 ' 2 ). 

The density condition allows to approximate rough solutions with smoother ones, 
locally uniformly in time, and is easily verified in concrete cases. 

Dispersive, smoothing and Strichartz estimates for a perturbed Dirac equation 
of the form 

iu t = Hqu + V{x)u 

were obtained earlier in [71 [HI [31 H] , for a general potential V = V* G C 4x4 satisfying 
suitable smallness and decay conditions. In those works we used a perturbative 
approach, relying heavily on spectral methods. Here we follow a different approach, 
based on multiplier methods, with two major advantages. First, we can partially 
overcome the smallness assumption; and second, the assumptions are more natural 
from the physical point of view since they are expressed in terms of the magnetic 
field 

B = curl A, 
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which is the physically relevant quantity. Actually, all assumptions are in terms of 
the quantities 

B T = ^AB, d r B = (d r B\d r B 2 1 d r B 3 ), 
\x\ 

which are, respectively, the tangential component and the radial derivative of the 
field B. Moreover, we establish for the first time a virial identity for the perturbed 
Dirac equation, which has several application not restricted to smoothing properties 
of the solution. Multiplier methods in relation with weak dispersion properties have 
a long story, starting from Morawetz [15] for the Klein-Gordon equation and [6], 
[TS] . [2U] for the Schrodingcr equation, and adapted to more general situations in 
[16j , [17] . Potential perturbations for the Schrodingcr equation were considered in 
PP, [2J, while the magnetic case was studied in [TT]. The perturbed Dirac equation 
was studied in [7], [5], [3] and [T] using spectral instead of multiplier methods, which 
are applied for the first time here. 

The paper is organized as follows. In the rest of the Introduction we shall de- 
scribe the main results of the paper, namely a general virial identity and optimal 
smoothing and Strichartz estimates for the Dirac equation perturbed with a mag- 
netostatic potential. Sections 2, 3 and 4 are devoted to the proofs. In the Appendix 
we prove a magnetic Hardy inequality (Proposition II. 5[) which is elementary but 
has maybe an independent interest. 

1.1. Virial identities. The Dirac operator does not have a definite sign and this is 
a substantial difficulty for a direct application of multiplier methods. To overcome 
it we shall resort to the squared Dirac equation 

(id t -H){id t +H) = (~d tt -H 2 ). 

Thus we are reduced to study a diagonal system of wave (Klein-Gordon) equations 
of the form 

(1.5) u tt (t,x)+Lu(t,x)=0, L=(m 2 -A)I 

with u = u(t, i) : I x R 3 -> C 4 . Our first result is a formal virial identity for solu- 
tions of a general system of wave equations like (|1.5p . with L being any selfadjoint 
operator on i 2 (R"; C k ). In the following, round brackets 

(F,G)= I F ■ Gdx 

denote the inner product in L 2 (R"; C k ), while [S, T] = ST- TS is the commutator 
of operators. 

Theorem 1.1 (Virial identity for the wave equation). Given a function <j> : K" — > 
M, define the quantity 

(1.6) 6(f) = (4>u tl ut) + 3? ({2<t>L - L<j>) u, u) . 
Then any solution u(t, x) of (|1.5[) satisfies the formal identities 

(1.7) e(t) = K([x,0Kut) 

(1.8) Q(t) = ~([L,[L,ct>}}u,u). 

We recall that the virial identity for a magnetic wave equation first appeared 
in [TTJ1 [TT] ; the abstract formulation given here can be applied to more general 
equations including magnetic Dirac systems. Thus, as an application of (|1.7p . (|1.8|) 
we obtain: 
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Theorem 1.2 (Virial identity for the Dirac equation). Assume the operator Ti 
defined by (|1.3|) satisfies assumption (A) and let <f> : K 3 — > R be a real valued 
function. Then any solution u(t,x) of (|1.4[) satisfies the formal virial identity, 
where DB = [9ji?i]i.j = i,3 and S = ja A a is the spin operator: 

(1.9) 

2/ V A uD 2 (jN~Xu~-)- I |u| 2 A 2 + 23 / u<f>'B r -Va~u + 2 [S ■ (DBWcj)) u] ■ u 

Remark 1.1. In the following, we will always consider radial multipliers 0, in which 
case the last term at the left-hand side of (|1.9p simplifies to 

3 ■ DBV<f> = <f>'S ■ d r B. 

In order to deduce a smoothing estimate from the virial identity, it will be 
necessary to impose suitable smallness conditions on the components of B appearing 
in (|1.9[) . for which no natural positivity assumption holds in general. 



1.2. Local smoothing for the magnetic Dirac Equation. As a first applica- 
tion of the virial identity, we prove some smoothing estimates for the magnetic 
Dirac equation. For / : M 3 — > C, denote by 



sup l/lllij = ( ( + (sup \f\fdr 

\x\—r \ J \x\— r 



\\J lli?L~(S r ) 

\x\— r \JQ \x\—r 

Moreover, by V^it and V^it we denote, respectively, the radial and tangential 
components of the covariant gradient = V — iA(t, x): 

(1.10) V>=^-V A u, S7 T A u = V A u-^V r A u, 

\x\ \x\ 

so that 

|V>| 2 + |V>| 2 = |V AU | 2 . 
We can now state our main result. 

Theorem 1.3. LetJi satisfy the self-adjointness assumption (A). Let B = curlj4 = 
B\ + B2 with B2 £ L°°(W l ) and inroduce the quantities 

Co = |||afBi|U~(R..), C\ = \\\^B T \\ LlL ~ {Sr) , C 2 = \\\x\ 2 d r B\\ LlLx(Sr) . 

We shall assume the smallness conditions 

1 



(1.11) C °<i' C 2 + 3C 2 + C 1V /C 2 + 6C 2 < 1 
and that the L°° part of B is absent in the massless case: 

(1.12) m = => B 2 = 0. 
TTien /or a// / G L 2 , the following estimate holds: 

(1.13) sup-/ / |e iw /| 2 aMi<||/||£ 2 . 

R>0 J-00 J\x\<R 
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Assume moreover that the second inequality in (jTTTTj) is strict; then for any f G 
D{TL) the following estimate is true: 

(1.14) 

sup - / / \V A e m f\ 2 dxdt + \\e m f\\ 2 LXL , 

— : ——dxdt + sup / / le^/l 2 ^ < ||W/||| 2 . 

Fl fl>0 J-oo J\x\=R 

Remark 1.2. Notice that all the assumptions in the smoothing Theorem 11.31 are 
expressed in terms of the magnetic field B, and consequently the gauge invariance 
of the result is preserved. 

Example 1.4. Explicit examples of magnetic fields satisfying assumption p. lip 
are of the following form 

(1.15) „(£)£ + eBW , 

where u is a smooth function on the unit sphere, while e is sufficiently small and 
B : K 3 -> M 3 satisfies 

1 .„ „, . 1 

for some S > 0. 



1^)1 ^ ^iTp+i ' \ drB \ ^ M 3-S 



In the proof of Theorem II. 3 1 we shall use the following Hardy inequality for the 
magnetic Dirac operator, proved in the Appendix at the end of the paper. Compare 
with [5] for parallel results. 

Proposition 1.5. Let B = curl A = B\ + B-2 and assume that 

(1.16) C = \\\x\ 2 B 1 \\ L ^ {K a) <oo, ||B 2 |U=o(R") < oo 

Then, for any f : R 3 — » C such that Tif € L 2 , and any e < 1, the following 
inequality holds: 

<L - ' I l/l 2 + - C) I ffi + « / IVx,P < (i + i^i / Wt 



4 

When to = 0, B-i = 0, the right hand side is to be intepreted simply as J \Hf\ 2 dx. 

1.3. Strichartz Estimates. As a natural application of the local smoothing esti- 
mates, we now derive from them the Strichartz estimates for the perturbed Dirac 
equation. We recall that the solution u(t, x) = e ltT> f of the free massless Dirac 
system with initial value u(0,x) = f(x) satisfies 

(1-18) \\e itv f\\ < \\fh*, 

LPH q q " 2 

for all wave admissible (p, q) 
2 2 2 

(1.19) - + - = -, 2<p<oo, oo>q>2 

p q 2 

while in the massive case m^Owe have 
(1-20) WefiV+nfifW i_i i < ||/ 



LPH. 



q p 2 



L- 



for all Schrodinger admissible (p, q) 
2 3 3 

(1.21) - + - = -, 2<p<oo, 6 > g > 2. 

p q 2 
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For a proof of these estimates see [5] • In the perturbed case we obtain exactly the 
same results: 

Theorem 1.6. Assume Ti and T>a satisfy (A). Moreover, assume that (|l.lip . 
(fTTT^T) hold and that 

(1.22) ^2 J sup \A\ < co. 

Then the massless perturbed flow satisfies the Strichartz estimates 
(1-23) lle^/H x , <\\f\\ L3 

L?H q q v 2 

for all wave admissible couple (p,q), (in particular, p ^ 2), while in the massive 
case we have, for all Schrodinger admissible couple (p,q), 

(1-24) \\e m f\\ i.i.x < ||/|| i2 (m^O). 

-j q p Z 



2. Proof of the Virial Identities 

Proof of Theorem \l.l\ The proof relies on a direct computation. By equation (|1.5p 
we have 

(2.1) ^( ( fm t ,u t ) = -2^((f>Lu,u t ), 
and 

(2.2) — 3?((2<£Z - L<j>)u, u) = ^St((<j)L + L<f>)u, u t ), 

since <^ and L are symmetric operators. Summing (|2.1[) and (|2.2[) we get (|1.7[) . An 
additional differentiation gives 

(2.3) 1 K([L ' Ut) = ^ Uu Ut) ~ 
Since [X, 0] is anti-symmetric, we have 

(2.4) U([L,<j>]ut,u t ) = 
and also 

(2.5) - SR([L, flu, Lu) = -i {([£, flu, Lu) + (Lu, [L, flu)} = ~([L, [L, fl]u, u). 
Identities |2~3]) . flU) and (HHJ) give JTHJ). □ 
Proof of Theorem II. £1 Let u be a solution to equation p.4p . Using the identity 

= {id t - H){id t + H)u = {-d tt - H 2 )u, 
we see that u solves a Cauchy problem of the form (|1.5|) : 

'u tt + H 2 u = 

(2.6) <u(0) = / L = H 2 . 

u t {<S)=iHf 

In order to apply (|1.7[) . (|1.8p . we need to compute explicitly the commutators 
appearing in the formulas with the choice L = Ti 2 . 

In the following we shall need the spin operator S, defined as the triplet of 
matrices 

_ i i , , 

b = -a A a = -{a 2 ct 3 - a 3 a 2 , a 3 ai — aia 3 , a±a 2 ~ a 2 ai)- 

We also recall the formula 

(2.7) (a-F)(a-G) = F ■ G + 2iS ■ (F A G) 
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which holds for any matrix-valued vector fields F = (F 1 , F 2 , F 3 ), G = (G 1 , G 2 , G 3 ), 
with F l ,G J G A'J4x4(C) (see [TO] for an extensive list of algebraic identities con- 
nected to Dirac operators). Thus expanding the square 7i 2 we have 

(2.8) H 2 = U\- H (a ■ A) - (a ■ A)H + (a • A) (a ■ A), 
where the unperturbed part is precisely 

(2.9) H 2 = (m 2 - A)/ 4 , 

and I4 denotes the identity matrix. Using (|2.7p we compute 
(2.10) 

-Ho{a-A)-(a-A)H + (a-A)(a-A) = i(V-A)+i(A-W) + \A\ 2 -2S-(V Ai + iAV). 
Now observe that 

(2.11) V A A + AA V = curLA = B 

(2.12) -A + i(V • A) + i(A • V) + |yl| 2 = (iV - A) 2 = -A A . 
In conclusion, by $n$, (|2TO|) . (j2~TU|) . (j2~TT|) . ([2"!^)) we obtain 

(2.13) H 2 = (m 2 - Aa)/ 4 -2S-B. 
Analogously, in the massless case we have 

(2.14) V 2 A = -A A h - 25 • B. 
Hence the commutator with reduces to 

(2.15) [H 2 , 4>] = [m 2 , 0] - [A a, 4>] - 2[S -B,<j>] = —[A a, </>]. 
Using the Leibnitz rule 

V A (fg) = gV A f + fVg, 

we arrive at the explicit formula 

(2.16) [H 2 , (/)} = —[A a, (/)} = -2V0 • V A - (A0). 
Recalling (|1.6p . I|1.7|) . wc obtain 

(2.17) e(t) = -$t(J u t (2V0- Vl^ + uA0)^ 

We now compute the second commutator between 7i 2 and 0. By (|2 . 13[) . (|2. 16[) 
we have 

(2.18) [U 2 ,[H\4>]] = [AaA^aA\] + 2[S-B,[A a ,4>]]. 
The term involving the magnetic Laplacian gives 

(2.19) (u, [A A , [A a , 4>]]u) =4 f V a uD 2 4>Va1I~ f \u\ 2 A 2 

+ 43 / u<j) f B T ■ Va~H 

(see formula (2.18) in [TT] with V = 0). By (|2~Td| . the last term in (|2~T5|) is equal 
to 

[5 ■ B, [A A , 4>}] = 2[S -B,V<f>- V A ] = 2(5 • BV(j) ■ V A - V<£ ■ V A S ■ B). 
Both <p and the components of the field B are scalars; moreover, we have 

[B,V A ]=-DB, 

where DB denotes the (matrix) gradient of the field B, and in conclusion 

(2.20) 2 [5 • B, [A A , <j>]] = -45 • (DBV(f>) . 

Finally, identity O follows from l[L"g|) . (j2~T7)) . (|2J8]) . ((2~T9j) and (j2~20l) . □ 
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3. The smoothing estimates 

The formal computations leading to the virial identity p.9[) make sense for suf- 
ficiently smooth solutions u £ C(R, if 3 / 2 ) and the choice of multiplier <j> we make 
below. Thanks to the density assumption (A), if we approximate data / E L 2 (resp. 
D(H)) with fj £ C^°, the corresponding solutions Uj = e ltn fj will converge to the 
solution u = e itn f in C([-T,T); I?) (resp. C([—T, T];D(H))) for all T > 0. 

We shall apply identity 1|1.9[) to the solution u = e ttn f of the problem 

(3.i) h; n ; 

with an appropriate multiplier function <f>. 

3.1. Choice of the multiplier. Writing r = |a;|, we define <f> as follows (see |llj ) 

pr— \x\ 

4>o(x) = / (^o(s)ds, 



00 = $>M = 





Jo 

where 

r < 1 

> r > 1, 

and M is a positive constant we will choose later. We have 

r < 1 
r > 1 

while the bilaplacian is given by 

A 2 0oW = -An6 x=0 - <5 N=1 . 
Moreover, for any R > we define 

so that by rescaling we have 
(3.2) 



R 



(3.3) 



(3.4) Afo(r) 




i±2M, r >R 



1 



(3.5) AVfl(r) = -4n5 x=0 - ^<* W = 



R- 



We notice that cf>' R , 4> R , A(/>r > and moreover 



(3.6) sup^(r)<M+i sup0&(r)<-^, A<t> R {r) < 1 + 2M . 
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3.2. Estimate of the RHS in (|1.9[) . Consider the expression 

/ u t (2V0 ■ V A u + uA0) = (tit, 2V0 ■ V A u + uA(f>) L 2 

appearing at the right hand side in (|1.9[) . Using the equation, we can replace Ut 
with 

Ut = —iHu = —imfiu — iD A u. 
By the selfadjointess of j3, it is easy to check that 

Si [— im(0u, 2V(f> ■ V A u) - im(0u, Acjm)] = 0, 

so that 

$t(u t , 2V0 • V A u + uA0) = 2$s(D A u, V0 ■ V^u) + ^(V A u, A0u) 
and by Young we obtain 

|2 i IIWJ, V7 ,.l|2 i 1 IL, A J.II2 



(3.7) 



3?(/ ut(2V0- V a m + uA0) 



<|lpA«||i 3 +||V0-VA«||i a +i||«A 



Recalling (|3.6p . and using Proposition 11.51 with the choice e = 1 — 4Co, which is 
positive in virtue of the assumption Co < 4 , we have 



(3.8) 



||V</>- Vau\\ 2 L 2 < 



The third term in (|3.7p , can be estimated using 
([TTTTjl : 

4 



and again the Hardy inequality 



(3.9) 



< 



(l + 2M)\\V A u\\i 2 . 



1 - 4C 

Summing up, by Q3.7p . (|3.8p and (|3.9p we conclude that 



(3.10) 

for any t 6 



9t( I u t (2V0- V A u + uA0) 
\Jr 3 



< 



\\T>au\ 



L- ■ 



3.3. Estimate of the LHS in (|1.9p . In this subsection, we bound from below the 
time integral of the LHS in (|1.9p . 

We shall use the elementary identity 



(3.11) 



V A uD^<pW A u 



V>| 2 + <£"(r)|V>| : 



By ([3~TTj) . (pHj) . l[3"3]) and ([33]), for the first two terms at the LHS of (TS]) we have 
(3.12) 



1 



|2a2. 



\u\ 2 da(x), 



V A uD z (j) R V A u - - / \urA 



\V A u\ 2 dx + 2M [ ^f± dx + 2Au{tM 2 + ^ I 

iH J\x\<R \X\ ZR J\ X \ =R 

for any R > 0, where der denotes the surface measure on the sphere of radius R. 
For the perturbative term involving B T in (|1.9p , by the Holder inequality and 
we obtain 
(3.13) 



23 



t=T 



t=-T 



u(f)' R B T ■ V A udxdt 



1 



t=T 



> -(2M + 1) ( sup-^ / / \u\ 2 da(x)dt 

R>0 R Jt=-T J\x\=R , 



t=T 
= -T 



\V A u\ 2 , J 
, , dxdt 



\x\ 2 B T \\ L 2 L a 
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For the remaining term in Q1.9P , observe that the operator norm of the components 
of S= {S\S 2 ,S 3 ) is 



|<S*]|c 4 ->-c 4 — 



hence we can write 

rt=T 



(3.14) 2/ / \u\ 2 S- [V<f> R DB]dxdt 

Jt=-T JR 3 

>-3^M+0 ^supi^ ^J^Ju\ 2 da R {x)dt^j \\\x\ 2 d r B\\ L i L <~ {Sr) , 
since 4>R is radial. Now we introduce the norms 

\\u\\ 2 XT := sup f [ \u\ 2 dxdt 

R>0 H Jt=-T J\x\<R 

\\ u \\yt : = SU P 4o / / M 2 dcr(a;)cft. 

i?>0 -K Jt=-TJ\x\=R. 

Taking the supremum over R > in (|3.12[) and summing with (|3.13[) . (|3.14[) . we 
obtain 

(3.15) 2 [ [ V A uD 2 cj) R V A ~udt - \ [ [ \u\ 2 A 2 cf> R dt 

Jt=-T JR 3 2 J t= -T JR 3 

+ 23/ f u<f> R B T -V^ + 2 [ [ \u\ 2 S -[V<j) R DB}dt 

Jt=-T JR 3 Jt=-T JR 3 

> ||V^||| iT +Q-3^M + ij \\\ x \ 2 d r B\\ L ^ iSr ^j \\u\$ tT 
-(2M+1) N* *J\x\- k *V T A u\\ 2 L2 dtJ |||s|*S T || i , £ - (Sr )||«||y ) T 

/•t=T /-t=T 

+ 2M / |||a;|-'V>||i 3 dt + 27r / |w(i,0)| 2 dt. 

Jt=-T Jt=-T 

In order to deduce (|1.13|) . (| 1 . 14[> . we need to ensure the positivity of the right-hand 
side of (|3.15|) . Define p, q as 

P= ^£"jlN-'V>]|| 2 ^ , q=\\u\\ Y , 

while C\ , C2 are defined in the statement of the Theorem. Then we are led to study 
the inequality 

(3.16) 2Mp 2 + ( i - 3 ( M + i J C 2 ] g 2 - {2AI + l)C lP q > 



and it is immediate to check that (|3.16p holds for allp, q > and M — C\j (2\J C\ + 6C2), 
provided C\ and C2 satisfy p. lip . Thus, dropping the corresponding nonnegative 
terms, we arrive at the estimate 

(3.17) 2 f f \7 A uD 2 R V A Hdt - \ [ f \u\ 2 A 2 <f> R dt 

Jt=-TJR 3 ^Jt=~TJR 3 
r t=T r rt=T 



23 f [ u(f>' R B T ■ V A u + 2 [ [ \u\ 2 S -[V(j) R DB]dt 

Jt=-T JR 3 Jt=-T JR 3 



>^Au\\' XiT >-\\V AU \\' XiT 
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where in the last step we used the pointwise inequality |2?aw| < |Vau|. We now 
integrate in time the virial identity on [— T, T], and using (|3.17[) and (|3.10p we get 

(3.18) f T \\V A u\\ 2 x dt < \\V A u(T)\\h + ||Z>au(-T)||£,. 

J-T 

Notice that all the above computations do not depend on the sign of the mass m. 

Let us consider the range of D A , form Proposition ll.51 we have that for Cq < 1 /4, 
^ ker(D A ) so Ra,n(D A ) is cither L 2 if is not in the essential spectrum of D A or it 
is dense in L 2 . Now fix an arbitrary g £ R&n(D A ), there exists / € D(D A ) = D(H) 
with T> A f = <?; we consider then the solution u(t,x) of the problem 

iu t = —m(3u + T> A u, it(0, x) = f 

with opposite mass, and notice that u satisfies estimate (|3.18[) . If we apply to this 
equation the operator Da we obtain, by the anticommutation rules, 

i(V A u) t = f3m{V A u) + V A {V A u), (D a u(0, x)) = V A f 

or, in other words, the function v = T>au solves the problem 

iv t = Hv, v(0,x)=g => v=-e ltH g. 

Hence p,18p can be written 

(3-19) f \\vf x dt < \\v(T)\\ 2 L2 + \\v(-T)f L2 = 2\\g\\v> 

J-T 

and letting T — > oo we conclude that 



/oo 
\\e Un 9\\ 2 x dt<\\g\\ L2 
-oo 



which is exactly p,13[) for g G Ran(I?A), which is dense in L 2 . So density arguments 
provide ()1.13(1 . 

In order to prove (| 1 . 14[) . let us come back to (|3.15p . If wc take the strict inequality 
in (|3.16p . which is equivalent to assume the strict inequality in p. lip , for instance 
the minorant being some positiv e, we have that 

(3.20) 

t=T i- ^ p t=T 



V AuD' z (/) R y A udt - - I I \u\ 2 k?cj) R dt 

-T JR a 2 J t= -T JR 3 

+ 2s/ / u<f>' R B T -W A 1 u + 2 [ I \u\ 2 S ■[V<f> R DB]dt 

Jt=-T JR 3 ' Jt=-T JR 3 

2 r t=T r t=T 

>^Au\\ 2 XtT + e\\u\\ 2 YT + e I |||a;|-iv>]|| 2 ^ + 27r / \u(t,0)\ 2 dt. 



t=—T Jt=-T 



By this inequality and (|3 . 1 0[) we obtain 



fliv^lli + 4Ay + 4\xr^>\\b + 2ttK*, o)| s 



dt 



<\\V A u(T)\\l, + \\V A u(-T)\\l 2 . 
The right hand side can be estimated using the obvious inequality 

\\v A f\\h < \\Hf\\h, 

and the conservation of ||Wu(i)||j,3. In order to complete the proof of (|1.14p . it 
only remains to remark that the term Hw-IIl^l^ m t ne inequality is obtained by the 
term jf T \u(t, 0)\ 2 dt by translating in space the multiplier <fi. Letting T — > 00 wc 
conclude the proof. 
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4. Proof of Strichartz estimates 



We pass to the proof of Theorem 11.61 for the massless case H = T> A ■ We rewrite 
u = e ltVA f using the Duhamcl formula: 

(4.1) u{t) = e itv f + f e'C-^a • Au(s)ds. 

Jo 

The term e'* 15 / is estimated directly via (|1.18|) . For the Duhamel term, we follow 
the Keel-Tao method (see [Tl], [12]): by the Christ-Kiselev Lemma in [5], it is 
sufficient to estimate the untruncated integral 



At-*)v a . Au(s)ds = e 



HP 



-isP 



a ■ Au(s)d 



since we are only interested in the non-endpoint case. Again by p,18[) we obtain 



(4.2) 



itP 



e~ lsU a ■ Au(s)ds 



LpH, 



1 _1_ 1 
q p ^ 



< 



-isP 



a ■ Au(s)ds 



L- 



Now we use the dual form of the smoothing estimate (|1.13[) . i.e. 



(4.3) 



e~ lsU a ■ Au(s)ds 



< 



U T,2 L 2 



L^L 2 (\x\^23)- 



j 61 



Hence, by Holder inequality, assumption (|1.22p and estimate (| 1 . 1 3|) we continue the 

estimate as follows 

(4.4) 

and this concludes the proof of (|1.23|) . 

We pass now to the proof of (|1.24[) in the massive case. By mixing free Strichartz 
estimates with the dual of (|1.13|) . for the Duhamel term we obtain 



(4.5) 



AtH B 



Au(s)di 



for any Schrodingcr admissible couple (p, q), with p > 2. The endpoint here can be 
recovered by using exactly the same technique as in [T3], Lemma 3. The rest of the 
proof is completely analogous to the massless case. 

Appendix A. Magnetic Hardy Inequality for Dirac 

We now prove Proposition 1 1.51 Denote by (•, •) the inner product in L 2 (R 3 ,C 4 ), 
|| • || the associated norm, and observe that, due to the formula (|2.7[) . we have the 
relation 

\\V A f\\ 2 = (a ■ V A f, a ■ V A f) = -((a ■ V A ){a ■ V A )f, f) 

= -(V A f,f)-2i(S-(V A AV A )f,f), 

where S = |o A a is the spin operator. Writing for brevity = dj — iA^ , we can 
compute explicitly 

Viav a = {{d£, d£], [dt d£], [dt d£]) = iB, 

where B = curL4. Hence by the previous relation we obtain 
0< \\V A f\\ 2 = \\V A f\\ 2 + 2(S-BfJ). 
Notice that S is a triple of matrices of norm < 1/2, hence we can write 

|2(S.5/,/)|<C |||*r7l| 2 
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and this implies 

(A.i) \\v A f\\ 2 > IIV.4/H 2 - Coll^H 2 - ||b 2 ||l~||/|| 2 . 

Now we recall the magnetic Hardy inequality 

(A.2) \ J j^pdi < J \V A f\ 2 dx, 

which is proved in [TT]. To complete the proof of (|1.17|) . it is sufficient to write 



\\nf\\l, = (H 2 fj) = 

and use the preceding estimates. 
The e inequality is obtained using 



\\V A f\\ 



(l-e)i J^dx + e J \V A f\ 2 dx< J \V A f\ 2 dx. 
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